Spin squeezing and entanglement in spinor-1 condensates 
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We analyze quantum correlation properties of a spinor-1 (/ = 1) Bose Einstein condensate using 
the Gell-Mann realization of SU(3) symmetry. We show that previously discussed phenomena of 
condensate fragmentation and spin-mixing can be explained in terms of the hypercharge symmetry. 
The ground state of a spinor-1 condensate is found to be fragmented for ferromagnetic interactions. 
The notion of two bosonic mode squeezing is generalized to the two spin (17- V) squeezing within the 
SU(3) formalism. Spin squeezing in the isospin subspace (T) is found and numerically investigated. 
We also provide new results for the stationary states of spinor-1 condensates. 



I. INTRODUCTION 

The availability of atomic Bose-Einstein condensates 
(BEC) with spin degrees of freedom have stimulated 
much recent interest in their applications to quantum 
information physics. Atomic spinor-1 condensates (hy- 
perfine spin / = 1 for each atom) were first realized by 
transferring spin polarized BEC prepared in a magnetic 
trap into a far-off resonant optical trap [3] . More recently, 
the all optical route to BEC has opened the door for their 
wide accessibility In a simple treatment, a spinor-1 
condensate can be described by a three-component order 
parameter, one for each of the Zeeman component of the 
hyperfine manifold. Early theoretical studies have clar- 
ified rotationally invariant descriptions including elastic 
s-wave collisions || |], pj, r| m] . Recent investigations re- 
veal that such a system also possesses complex ground 
state structures and can exhibit novel dynamical effects 
J?], H, ^| , such as fragmentation ]ic| , spin-mixing and en- 
tanglement d 0, y| . 

In this paper, we further explore quantum correlation 
properties such as spin squeezing and entanglement of a 
spinor-1 condensate [n], 121. For a spin half (/ = 1/2) 
atomic system, a rotationally invariant Hamiltonian is 
known to not induce spin squeezing as the total spin is 
conserved jl3). For a spinor-1 condensate, it was found 
that its Hamiltonian becomes rotationally invariant if the 
single (spatial) mode approximation is made to its order 
parameters, i.e. assuming jj>m,f (^) = 4 > {^) a m f with the 
same mode function <j>(r) |7]]. a mf is the annihilation 
operator for atoms in Zeeman state mj = ±, 0. However, 
as we show in this work, various nonlinear processes do 
occur within different subspaces of the full SU(3) struc- 
ture of a spinor-1 condensate, e.g. we find the existence 
of spin squeezing in the isospin sub-group jl3 . 

The generation and characterization of squeezing and 
entanglement of Bose-condensed atoms has recently 
emerged as an active research area. Historically, atomic 
squeezed states were first considered for a system of two 
level atoms. Even in this simple SU(2) case, it was found 
that some operational definitions of spin squeezing can 
become system dependent. Depending on the context in 
which the concept of squeezing is applied, different defi- 
nitions arise uM . Squeezing in atomic variables was first 



introduced through reduced fluctuations in the atomic 
(Pauli) operators of the system |l5| such that atomic res- 
onance fluorescence in the far-field zone is squeezed. In 
this case, it is useful and convenient to define the atomic 
squeezing parameter according to, 



£ h = AJ i /^\(J j )/2\, i^je(x,y,z). 

This definition can be essentially read off from 
Heisenberg uncertainty relation AJ^AJj > |(Jfc)/2| 



(1) 

the 
for 

the collective angular momentum components of the two 
level atomic system. When in such a squeezed state 
£/j < 1, the quantum fluctuation of one collective angular 
momentum component becomes lower than the Heisen- 
berg limited value at the cost of increased fluctuation in 
the other component. The general family of two level 
atomic states satisfying this criterion was found to be 
Bloch states, or SU(2) coherent states. These "squeezed 
states" are obtained by simply rotating in space the col- 
lective (Dicke) state |J, ±J) |l(|. Soon it became clear 
that neither £/, nor the Bloch states are that useful in 
other applications of atomic squeezing. In particular, for 
Ramsey oscillatory field spectroscopy, a new squeezing 



parameter 17 



& = V2JAJ±/\(J)\ 



(2) 



is called for with Jj_ the angular momentum component 
normal to the (J), i.e. in the direction of the unit vec- 
tor n along which A(n • J) is minimized. The squeezing 
condition £r < 1 is not straight forwardly determined by 
the Heisenberg uncertainty relation. Instead, it is defined 
by requiring the improvement of signal to noise ratio in 
a typical Ramsey spectroscopy. It was later shown that 
the same criterion is also applicable for improving the 
phase sensitivity of a Mach-Zehnder interferometer juj . 
£r = 1 for Bloch states or SU(2) coherent states. An 
independent refinement of £h was suggested by Kitagawa 
and Ueda [ fl3| to make it independent of angular momen- 
tum coordinate system or specific measurement schemes. 
They emphasized that collective spin squeezing should 
reflect quantum correlations between individual atomic 
spins and defined a squeezing parameter 



e 9 = aj±a 



(3) 
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to measure such correlations. The factor J/2 in the 
denominator represents the variance of a Bloch state, 
which comes from simply adding up the variance of each 
individual spin (1/2). When quantum correlations ex- 
ist among different atomic spins, the variance of certain 
component of the collective spin can become lower than 
J/2. This leads naturally to the criterion £ q < 1 for spin 
squeezing. We note that this definition is directly re- 
lated to the spectroscopic definition as £r = ( J/\( J z )\)£ q - 
More recently, a particular type of quantum correlation, 
namely the multi-particle entanglement, becomes impor- 
tant for quantum information physics. A more stringent 
criterion for atomic squeezing which combines the quan- 
tum correlation definition with the inseparability require- 
ment of system density matrix is given by S0rensen et al. 



(2J)A(ni • Jf 



< 1, 



(4) 



(rf 2 • J) 2 + (r? 3 • J) 2 

with the n, being mutually orthogonal unit vectors. £ e is 
in fact identical to £r along the direction (ni • J) = 0. It 
was proven rigorously that when £ e < 1, the total state 
of the N two level atoms becomes inseparable, i.e. en- 
tangled in a general sense. All three definitions above 
applies to a two component (two level) atomic system. 
Many complications arise when attempt is made to ex- 
tend spin squeezing to a spinor-1 SU(3) system. Under 
certain restrictive conditions, £ e has been used recently 
to discuss a two mode entanglement in a spinor-1 con- 
densate 0. 

A related problem to spin squeezing is its efficient gen- 
eration and detection. In accordance with their respec- 
tive definitions for SU(2) systems, several physical mech- 
anisms have been proposed along this direction. Kite- 
gawa and Ueda considered a model Hamiltonian Hku — 
TiX J 2 that can be realized via the Coulomb interaction 
between electrons in the two arms of an interferometer 
fl3|| . Barnett and Dupertuis suggested that spin squeez- 
ing can be achieved in a two-atom system described by 
Hbd = ih{g* J1+J2+ — h.c.) [[18). Use of a pseudo-spin 
two component atomic condensate system has also been 
suggested |1| . Recently, two different groups considered 
atomic (spin) squeezing and entanglement in a spinor-1 
condensate, under the assumption that one of the com- 
ponent is highly populated such that quantum properties 
are important only am ong the remaining two sparsely 
populated components Jll| , [12]] . Our aim in this study is 
to remove such a restrictive condition, and consider the 
full quantum correlations within a spinor-1 condensate. 

Let us consider a general three component system la- 
beled by i,j G {+,—,0}. For spectroscopic and inter- 
ferometric applications the observables of interest are 
the relative number of particles Ni — Nj (particle par- 
titioning) and the corresponding phase differences <pi — 
with their measurements limited by noises SNij — 
([A(Ni - Nj)} 2 } 1 / 2 and S&j = ([A(& - ^)} 2 ) 1/2 - For 
a two-component system, the particle partitioning be- 
comes the collective angular momentum projection as 



N + — AL. = 2 J z in the standard Schwinger representa- 
tion; the relative phase becomes the corresponding az- 
imuthal phase <f> z = <p + — <j>_, which is conjugate to 
J z . Quantum mechanically they satisfy [J z ,<f> z ] = i. 
Thus from 5J z S<j) z > 1 and 5N+_ = 2(AJ 2 ) 1 / 2 , we 
find 5(f> + - w (AJ 2 y/ 2 /\(J x )\ g§. Therefore, for spin 
squeezed states, one achieves higher angular resolution 
and reduced particle partitioning noise. In a three com- 
ponent system, we can similarly associate the three num- 
ber difference Ni — Nj with three subspace pseudo-spins 

(each of spin-1/2) U, f, and V such that N + - N- 
■11'.. N+- N = 2U 3 , N_-N = 2U 3 . The phase differ- 
ences can then be similarly expressed in terms of com- 
ponents U X y, V Xj y, and T Xj y. When demanding noise 
reduction in such a SU(3) system, we need to consider 
squeezing in the three spin-1/2 subsystems. One may 
naively expect that results from the above discussed 
SU(2) squeezing can be applied to each of the three sub- 
systems, and collectively, one can simply demand that 
£ e < 1 to be satisfied simultaneously. In reality this 
does not work as the three spin-1/2 subsystems do not 
commute with each other. This is also the fundamen- 
tal reason that makes it difficult to generate and de- 
tect quantum correlations in a full SU(3) system. Fur- 
thermore, due to the above non-commuting nature, the 
three SU(2) sub-spins cannot be squeezed independently 
of each other. Previous discussions of a spinor-1 conden- 
sate entanglement are always limited to just one SU(2) 
subspace, usually in the limit Nq ~ N, i.e. one mode 
is highly populated. Approximately, this limit destroys 
the underlying non-commutative algebra among (U, V, T) 
and simplifies the problem to that of a usual two-mode 
SU(2) system. 

One of the major results of this paper is that the ef- 
fective Hamiltonian of a general spinor-1 condensate can 
be decomposed as 



H 



^XKijTl + %Xbd{U + V + + h.c), 



(5) 



which involves both the Kitagawa-Ueda (KU) and 
Barnett-Dupertuis (BD) type of spin squeezing simul- 
taneously. In other words, all three fictitious spins can 
indeed be found squeezed in a spinor-1 condensate, as 
the above two distinct nonlinearities commute with each 
other, and therefore squeeze all three SU(2) subspaces 
simultaneously. We find that the BD type interaction 
dominates when A'o is large; while in the opposite limit 
the KU type squeezing governs. For intermediate values 
of A'o it is necessary to consider a generalized spin squeez- 
ing for the three mode spinor-1 system. To achieve this, 
we provide a new criterion for the U-V two spin squeez- 
ing based on reduced quantum fluctuations imposed by 
the BD type nonlinearity. When such a condition is satis- 
fied, the state of a spinor-1 condensate as a macroscopic- 
coherent quantum object becomes useful for three-mode 
spectroscopic and interferometric applications. We fur- 
ther show that this condition also corresponds to a two- 
mode entanglement in terms of the Holstein-Primakopf 
bosonic modes, and it reduces to previous results in the 
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large Nq limit |Tl], [TJ| . Squeezing in T spin is particularly 
useful for quantum information applications based on col- 
lective (Dicke) states | J, J z ). These states are in fact sta- 
tionary in a spinor-1 condensate and can be manipulated 
via external control fields 0. Since J z = N + — N- = 2T 3 , 
such T-squeezed states ensure well-defined Dicke states. 

In our study of spin squeezing in a spinor-1 conden- 
sate as outlined in this paper, we present a systematic 
approach by recognizing the (U, V, T) pseudo-spin sub- 
spaces as the Gell-Mann (quark) realization of the SU(3) 
algebra |20[ . Similar recognitions are found useful in 
the recent discussions of quantum and semi-classical dy- 
namics of three coupled atomic condensates pH , where 
the BD-type two-spin squeezing nonlinearity was absent. 
Earlier investigations of three level atomic systems also 
made efficient use of the density matrix and expressed 
atomic Hamiltonians in terms SU(3) generators (22[ Eg) . 
The main difference between our approach (on spinor-1 
BEC) and those earlier studies are the enveloping Weyl- 
Heisenberg algebra of the bosonic operators, which leads 
to subsequently much larger Hilbert space of the sys- 
tem. In addition to spin squeezing, we also investigate 
other quantum correlation effects, e.g. condensate frag- 
mentation with the new theoretical framework. We show 
that previous theories based upon the SO (2) rotational 
symmetry group cannot give a decomposition of angular 
momentum operator with nonlinearities that could eas- 
ily be considered for spin squeezing. For instance, neither 
H oc L 2 -2N § nor H oc N(N—1) — A^A [j§ can lead to 
any simple recognition of the nonlinear coupling among 
various spin-components. 



II. SU(3) FORMULATION FOR SPINOR-1 BEC 

Under the single-mode approximation [Q, a spinor-1 
condensate is described by the Hamiltonian 



Y+2a 



H = fiN — X' S N{N + 1) + \' a {L 2 - 27V), 



(6) 



where N = n + + n_ + no is the total number of atoms 
and the collective angular momentum (L) with the famil- 
iar raising (lowering) operator L + = \/2(a+ao + c4 a -) 

(L_ = L z = n+ — n_, and rii — at a*. u is 

the chemical potential. a are renormalized interaction 
coefficients, they are related to various s-wave scatter- 
ing lengths and <p(r) || fjj. The validity of the single- 
mode approximation is now reasonably well understood 
jj], [y], and for ferromagnetic interactions, it is in 
fact exact as shown recently [24}| , a±.o can form a sim- 
ilar Schwinger representation of SU(3) in the following 
manner, 



T+ = 



a+a_, r 3 = -(n + -n_), 



a±a , 



v 3 \_i 

U 3 J 2 



N 



n + + n- + no, 



(7) 

(8) 
(9) 




FIG. 1: The action of Gy in T3 — Y space. Any point is cou- 
pled only to its next nearest neighbors along T3-axis through 
a two step process V+U+ on the V-line with T3 unchanged 
in the end. Note that V+ and U+ commutes with each other 
and the conjugate process is also shown. 



Y = -(n+ + n_ - 2n ). 



(10) 



The linear combinations X± ± X T (X = T, U, V) to- 
gether with T3 and Y resemble the set of eight gen- 
erators of SU(3) in spherical representation. T±^,U±^ 
and V± y s fulfill commutation relations [A + ,A_] = 2A3 
and [X 3 , X±] = ±X± of the SU(2) algebra. We call T- 
operators the isospin and Y operators the hypercharge 
only because of their formal resemblance J^TJ. U and V 
subalgebras will be called U- and U-spin, respectively. 
We then have L + = V2(V + + U-), L z = 2T 3 , and 



L 2 =4T 3 2 + i(/V-e + )(/V- e _) 



-2(Y -Yof + G Y , 
= 2{V+U+ + h.c), 



(11) 
(12) 



with e± = -3/2 ±V2, and Y = -N/6-1/4. For N > 1 
this gives Yq ~ —N/6, which corresponds to no — N/2. 
Using [T 3 ,V±] = ±V±/2 and [T 3 ,U±] = TU±/2 we find 
[T 3 ,Gy] = consistent with [H,L Z ] = 0. Hence the 
Hamiltonian (^|) separates into three commuting parts 
H = H N [N] + H T [T 3 ] + H Y [Y]. To our knowledge, this 
decomposition have not been discussed before. In Ref. 
pit , a model of H = x{T 3 + 3Y 2 ) has been considered 
for both the quantum and semi-classical dynamics of Y 
as well as for SU(3) coherent states. We note that the 
decomposition ( |l l| ) differs from the Casimir relation for 
the two mode case pSj ]. In fact, with the spin singlet pair 
operator A — {oq — 2a+a_)/v / 3 as defined by Koashi and 
Ueda H, we find L 2 = N(N + 1) - A* A. 

Denote the simultaneous eigen-states of commuting op- 
erators (N,Y,T 3 ) as |7V,T 3 ,y), we find 



N,T 3 ,Y) 




Y 



N 



±T 3 



Y 
~2 



\\\N,T 3 ±-,Y 
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i.e., Gy only couples next nearest neighbors along the 
Y-axis through off-axial hopping as depicted in Fig. 1. 
Perhaps it is not surprising that operators T±, U±, and 
V± are simply the off-diagonal elements of the single par- 
ticle density operator = a)^a Ul while N, T3, and Y 
are related to the diagonal elements. T±, V^, U± all raise 
and lower the T3 value by (1 or 1/2). 

As an example we consider the simple case of the X3 = 
block along the line in Fig. 1 of the Hamiltonian (||) for 
X' a < 0, the ferromagnetic case (asin 87 Rb §). The polar 
case of X' > (as in 23 Na BEC [Q) has been discussed 
in Ref. ||. Dropping the constant Hn[N] and write in 
units of |A„|, the Hamiltonian becomes 

H = 2(Y- Y ) 2 - 2t Y (\Y + 2)(Y\ + h.c.) , (13) 

where t Y = (N/3 - Y)(N/3 + Y/2 + 1). Following the 
tight-binding procedure for the restricted two-mode case 
(+ and — ) discussed in Ref. |J, its eigenstates can be 
found by determining the ip(Y) of \ip) — Y^y V'OOI-^O 
through a difference equation 

Eil>(Y) = 2(Y-Y ) 2 ^(Y) 

-2[t Y - 2 ^{Y -2)+t Y il>(Y + 2)]. (14) 

In the continuum limit and up to the first order in 
0(\Y\/N), the equivalent differential form becomes 



0.25 



E 



8Y 2 



1 U = -2 



OY 2 

Its ground state is therefore 



1 dtp 
YodY 



(Y-Ypf 
4Y 2 



i>. (15) 



E 

Y 



exp I 



V2(Y-Y f 



Y-Y„ 



(7r|Yb|/V2)V4 



\Y), (16) 



which gives a diagonal (p^ v ) with (no) = N/2 and 
(n±) = N/4, i.e. a fragmented state |9[ [To). To check 
the validity of this approximate analytical result, we also 
solved the same problem for N = 10 3 within the T3 = 
block by an exact diagonalization procedure. The results 
are compared in Fig. |^. We see that without any fitting 
parameter the analytical result agrees well with the exact 
numerical result. Since the Hamiltonian is block diago- 
nal in even hq and odd no spaces, we find two degenerate 
ground states with even and odd no components respec- 
tively. These even-odd ground states display opposite 
phases and can form a Schrodinger cat state |M . The ap- 
proximate result here applies for a value of even N, which 
leads to an even no within the T3 = block. Hence, only 
the even no block of the Hamiltonian is considered. We 
will also show below that under the single-mode approx- 
imation, the exact ground state for N 3> 1 is generally a 
fragmented state with (n ) = N/2. More detailed studies 
and implications of the cat-like ground state in the even- 
odd number blocks and with their respective phases will 
be explored further and results presented elsewhere. We 
note that such cat states separated in the angular mo- 
mentum L z have been found in studies of Josephson type 




600 



FIG. 2: The ground state expansion coefficients ij)(Y) = 
ip(no) as a function of no for TV = 1000 atoms in the T3 = 
block. The solid curve is the approximate analytical result Eq. 
(15) without any fitting parameters while the other curves are 
obtained by an exact diagonalization procedure. 



coupled condensates In a spinor-1 condensate 

as considered here, we hnd that the dynamical behavior 
can be characterized by Y — N/3 — no, which can be ex- 
pressed as Y = 2(1/3 + V3) /3. Since the azimuthal phases 
are conjugate to angular momentum z-projection opera- 
tors, the cat-like ground states predicted here resemble 
the angular momentum cat states of a two-component 
condensate in its conjugate phase spaces. Finally we note 
that the symmetry point Y"o can be adjusted by exter- 
nal control fields which contribute terms proportional to 
n±_o to the Hamiltonian, and can be absorbed into the 
(Y~Y ) 2 term through a new Y a as n± = N/3 + Y/2±T 3 
and n = N/3 - Y. 

We have now seen that the Hamiltonian of the system 
describes an effectively one dimensional dynamics along 
the Y-axis, similar to that of a diffusive random walk pro- 
cess but now with an attractor (for X' a < 0) Y). Hence, 
we expect Y) influence population dynamics in a similar 
manner it effects the fragmentation. For T3 = 0, it is 
known that populations oscillate around time-averaged 
values no = N/2 and n± = N/4 which are the same as 
the results we found for fragmented ground states. We 
conclude that steady-state values of population oscilla- 
tions as well as fragmentation is determined by the hy- 
percharge symmetry point Y), which can be shifted by 
external fields. 



III. TWO-SPIN AND ISOSPIN SQUEEZING 

The form of Gy suggests the existence of two mode 
squeezing as was also noted recently by Duan et al. [jll| , 
who studied a spinor-1 condensate initially prepared in 
the Fock state with only m,f = state populated. During 
the time when the total number of excitations into states 
rrif — ±1 are negligible, the spin mixing term (Gy) in 
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the Hamiltonian simply reduces to a two-mode squeezing 
nonlinearity via (no) (a + aL +h.c). This creates a contin- 
uous variable type entanglement, or mode-entanglement 
in the second quantization form. In order to relate it 
to measurable spectroscopic spin squeezing and particle 
entanglement, Ref. (T^j first showed that in the low ex- 
citation limit, the two-mode entanglement criterion can 
also be expressed in terms of spin squeezing parameters 
for L x<y . In order to use the two-level SU(2) definition 
for spin squeezing of L xy , new pseudo-spins J± was in- 
troduced within the the two level subsystems | + 1) ± | — 1) 
and |0). They found that when L z as 0, the system 
Hamiltonian becomes effectively H = X a (L 2 — 2 A) w 
K(L X + Ly) ~ (J\ x + Jl y ), which causes each spin 1/2 
subsystem to be squeezed via the single axis twisting 
scheme. For the independent single-axis twisting scheme 
to work efficiently in achieving substantial spin squeez- 
ing, the commutator [J+ x / y , J- x /y] — (T- — T+)/4 oc T y , 
needs to be small. Hence, squeezing in the isospin is es- 
sential to achieve this two-mode squeezing goal. Without 
it, large quantum fluctuations in T y would destroy the 
two mode squeezing. Unfortunately, both the relation- 
ship between the two-mode squeezing and spin squeezing 
as well as the interpretation in terms of a dual single 
axis twisting fails to be adequate for higher excitations 
under more realistic situations. Indeed, for the extreme 
opposite case of % < n± ~ A/2, the Hamiltonian de- 
scribes a single-mode amplitude squeezing as it reduces 
to Gy ~ (aj 2 + h.c). Anywhere in between of these 
two extreme limits, we propose a new type of squeezing, 
the two-spin squeezing as a generalization of single spin 
squeezing by taking into account quantum correlations 
for mode-entanglement applications. We first note that 
the two extreme types of squeezing in Gy can be han- 
dled at arbitrary levels of excitation by introducing a new 
two-spin squeezing operator via 



two-mode squeezing or amplitude squeezing in the ap- 
propriate no limits. At low excitations when no ~ A, 
we have A3 w —A/2, S x w (no), and A^ w 0. In 
this case, A_ « y/ (no)a x and Gy = 2(no)(a) v a\ l + h.c.) 
demonstrates the two-mode [SU(1,1)] squeezing as in Ref. 
pT[ . In the large Hq scheme of Ref. O] , such modes are 
sparsely populated since, a'J_a' + = n + (l + no)/no- In the 
opposite case of large n± , we are in the strong excitation 
regime with N v « 1,S V / U ~ (n±) which gives effective 

modes to be a±ao/ \Jn± with large occupations. 

In order to define two-spin squeezing introduced via 
the A-operators in a similar way to the two-mode bosonic 
squeezing, we introduce Hermitian quadrature operators 



Q- 



K 
Q+ 

fir/2 



(e M C/_ + h.c.)/V2, 
(e ia V- + h.c.)/V2, 

{xz + xs)/2, 

( X «+V2 -XS +n/2 )/2. 



(19) 
(20) 
(21) 

(22) 



From J± = V - U and J ±3 = (3Y/2 ± T x )/2 we find 
= n(a) ■ J± with n(a) = (cos a, sin a, 0). If U and V 
were uncorrelated, their respective quantum noises would 
contribute to that of J± additively. Existence of quan- 
tum correlations between the U- and U-spins would re- 
duce the quantum fluctuation in J±. Thus, the (J±) spin 
squeezing is achieved by two-spin (U-V) squeezing. From 
[U-,V-] = 0, we find that 



(AQl) 2 + (AQ^ /2 ) 2 
= £ [(AX:) 2 + (AA»+"/ 2 ) 2 l 



with the U-V correlation function 



not 

uv 



-2ia 



(V + ,U+)/2 + c.c. 



(23) 



(24) 



A' 4 



V + U 4 



a\a_ 



n » n±, 
n < n± , 



K- = K\ 



(17) 
(18) 



with [K-,K+] — 2 A3. The squeezing mechanism in 
Hamiltonian (||) is now understood to be a general- 
ized Barnett-Dupertuis (BD) type squeezing via the 
V+U+ + h.c. nonlinearity in Gy. This is significantly 
more complicated than the two bosonic mode squeezing 
as the two spins U and V have a non-commuting al- 
gebra. The mode-entanglement of approximate bosonic 
modes a' ± = a±a/ J yj (no) of Ref. |ll| can in fact be gen- 
eralized to mode-entanglement between exactly bosonic 
Holstein-Primakoff modes |E3 a x — UlV defined through 
A+ = al^S x - N x and A 3 = N x -S x /2, in the spin S x /2 
realization of corresponding SU(2) algebras of U and V 
spins with N x — a x a x . The squeezing treatment with 
the exact bosonic modes ajj and ay remains to be more 
complicated than the usual two bosonic mode squeezing 
as it also suffers from the underline non-commutating al- 
gebra. This representation reduces to the usual SU(1,1) 



denote the correlations among U-V spins to the quadra- 
ture noise, which reduces the uncertainty bound when 
two spin squeezing occurs. We find a lower bound for 

the quadrature noise £ 0=1ti1) [(AX£) 2 + (AA a Q+7r/2 ) 2 ], 

by noting that [X", X« +7r/2 } = 2iV 3 , [A«,A^ + " /2 ] = 
-2iU 3 , and \(U 3 )\ + \(V 3 )\ > \(U 3 + V 3 ) \ = 3\(Y)\/2. We 
finally find 



(AQ°) 2 + (AQ^ + * /2 ) 2 > 3|<r)|/4 + C^ 



(25) 



Therefore, taking into consideration the important spin- 
spin correlation between different particles similar to the 
1/2 case OL we can introduce the U-V squeezing con- 
dition as 



So 



A(Q«) 2 



(AQ_; +7r/2 ) 2 



\(Y)\ 



<3/4, 



(26) 



similar to the continuous variable system |3C|| . This is the 
major result of our paper on the two [SU(2)] spin squeez- 
ing within the SU(3) of a spinor-1 condensate. The signif- 
icance of spin-spin correlation function to spin squeezing 
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and entanglement for a two-mode system was previously 
discussed in Ref. [[n] , where they showed that a negative, 
finite correlation parameter causes spin squeezing and 
entanglement of the atomic states. With the Holstein- 
Primakoff relations, it is straight forward to show this 
condition contracts into (£") 2 + (£" +7r / 2 ) 2 < 2 when 
no — > N. Thus Eq. (p6|) generalizes the two- mode entan- 



glement criterion (£") 2 + (£" +7r ^ 2 ) 2 < 2 at low excitations 
p| to arbitrary levels of excitation for two spin squeez- 
ing. For completeness, we note the squeezing parameter 
for J± -spins are O] 



N((AQ%) 2 



(Q± + * /2 ) 2 + (J±s) 2 '' 



(27) 



we have now greater freedom to consider a suitably 
prepared spinor-1 condensate to achieve many-particle 
and/or mode entanglement for quantum information ap- 
plications as well as various type spin squeezing for atom 
interferometry and spectroscopy applications in the long 
time limit with more macroscopic populations in all / = 1 
three component states can occur. In the limiting case 
discussed before either uq ~ N or n± ~ N is required 
to be large, the quantum states (modes) of interest are 
always sparsely populated. More generally, one can use 
Raman coupled laser pulses on a spinor-1 condensate to 
generate states with arbitrary populations in each mode 
and with arbitrary initial phases. This allows then for 
the consideration of stationary states in the fully quan- 
tum mechanical framework for their use in squeezing- 
entanglement applications. 



while the Heisenberg uncertainty relation gives 
AQJAQ± +7r/2 > |(J± 3 )|/2. For many particle en- 
tanglement of three level atoms, the criterion is given by 
either £± < 1. 

The U-V squeezing discussed above displays existence 
of nonlinear interactions within/among T, U, and V sub- 
space of (|). One may also contemplate for a one-axis 

isospin twisting (through T3) of the particular form of 

L 2 (In]). However, the dynamics of spinor-1 BEC be- a-di + a + a\) N \0, 0, 0)/VM, can be prepared with Ra 



man pulses [Q , where |0, 0, 0) is the vacuum in the Fock 
basis j no, n_, n + ) and otj — \aj\e zSj complex. Using 
to = n + — ri- , we can write 



IV. RESULTS AND DISCUSSIONS 

We now present some results on the numerical inves- 
tigation of isospin squeezing. If the condensate atom 
number N is fixed, a generic state |V>(0)) = (aoa + 



comes considerably more complicated because of off- axis 
hopping processes along the hypercharge axis (as in Fig. 
1). Due to the non-commutativity of sub-spin systems 
{U, V, T), squeezing and entanglement appears even with- 
out essentially any axis- twisting. In fact, even when 
T3 = 0, squeezing within the isospin subgroup can still 
happen as the U-V two-spin squeezing interaction would 
redistributes the noise also for the isospin subspace, in 
addition to the U-V spin space. To appreciate this fact, 
let us consider the rotation operator involving only U- 
V spins and employ the SU(2) disentangling theorem to 
obtain 



1^(0)) = y^^Nmk(a) 



2k, 



N — TO 



N + m 



where a — (ota, a-, a+), k = 0, 1, • • • , (N — |to|)/2 for 
even N + m, 2k = 1, 3, • • • , (N— |to|) for odd N + m, and 



N-m Af-r 

V'jvmfe = \/ C 2 N k C N * 2k afa_ 2 



R[C 



e^+(l + |7 ? | 2 )^e-"* i -, (28) where C™ = 



denotes the binomial coefficient. 



with 77 = CtanC/|CI- Using [V+, U-] = T+ and [V+,T + ] = 
[U-,T+] = 0, we find 



e r,L + _ e riV2V +e riV2U- e -r,T+/V2^ 



The basis transformation coefficients between angular 
momentum and Fock sates are available from Ref. 
written in more compact forms as 



(29) 



Hence, we arrive at 



\lm) =J2 G 



link 



, N — m , N + m 
2k,— k,— k ) (31) 



R[Q = e 



y/2viV +p V2riU- 



R T [rj\e 



-V2r)*U +e -V2ri*V- 



, (30) 



with a rotation operator within isospin space via Rt = 
exp {-T]T + /V2)(1 + h| 2 ) 2T3 exp {-r,*T^/^/2). This re- 
sult reflects the nature of Euler-angle rotations in three 
dimensions for a spin-1 system. We thus conclude that 
squeezing in J± through redistributing the noise via ro- 
tations is always accompanied by a redistribution of the 
noise in the isospin subspace. Squeezing and many parti- 
cle entanglement via the isospin can be checked using 
the usual spin squeezing criterion, which for both T- 
squeezing and the above derived U-V squeezing are in- 
dependent of their respective initial conditions. Hence, 



with 



Glmk — 2 fc s/ ^ 



(-ir 



Nlm 
kr 



and the symbolic notation 



Nlm 
kr 



xC, 



(l-m)/2-k+r 
l-2k+2r 



(N-m)/2-k ( jl- 



N-2k 



21 



(32) 



We note I = N, N - 2, • • • , N - 2[N/2] with [n] 
n,n — 1/2 for n = even, odd, and r = max[0, k 
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(I - \m\)/2],---,mm[k,(N - l)/2], m = 0, ±1, • • • , ±Z. 
k = 0,1, • • • , (N - \m\)/2 for / + m = even and 2k = 
1, 3, • • • , (N — | to |) for I + m = odd. The normalization 
is given by 

x (JV-0/2 

^2 = JJ C 2j C (N+l)/2-j- ( 33 ) 

Simpler analytic results exist for special cases, e.g. 

(j-Nmk — ^ V N ^N-2k l\l L '2N > l d4 i 

which takes an asymptotic form Gjvofc — ^/ Cff /2 N ~ 1 
when AT 3> 1. We also find 

AT/2 

K> = ^rE^ fe ( 2fc ) = iV A (35) 

fe=0 

the same result as obtained previously from (|l6|). 

These expressions allow us to express the initial state 
as |V>(0)) = E; m ^m(O)l^m), with 

V>zm(0) = Ajv m (a_,a + )E77 fc (d)Gjv mfe G imfc , (36) 

k 

a / ' n N-m (JV-m)/2 N+m/2 

Antu = y C 2JV a_ a + ! (37) 

and 77 = aQ/(2a_a + ). When 77 = 1, it becomes an eigen- 
state of Hamiltonian (^) as J2k GNmkGimk —> §ni- This 
generalizes the stationary state of J?J into the quantum 
regime. The condition for stationarity becomes 2Sq = 
5-+5 + and + = 1 (since |a | 2 + |a_| 2 + |a + | 2 = 
1). For the special case of Sj = and a_ = a + we obtain 
Uq = 1/2, in complete agreement with earlier results JtJ. 
By defining Pj = \aj\ 2 as spin component populations, 
we find that stationary states require Pq = 1/2 whenever 
P_ = P + . This is, however, not sufficient without estab- 
lishing the phase constraint found above, which becomes 
particularly useful as it provides for more freedom in state 
preparation using Raman coupled laser fields. As an ex- 
ample, we now consider isospin squeezing with the same 
form of initial states as in Ref. § for a = ^Po~e i9 / 2 
and a± = vT^Tb- This gives (T x (0)) = N(l - P )/2 
as the only non- vanishing isospin component (at t = 0). 
The population in the nif = component then acts as a 
knob between the two extreme squeezing type discussed 
earlier as well as between the Gy and T 3 terms. In the 
special case of within the T3 = block, we find that 
the dynamics of the system is determined only along the 
hypercharge Y axis. Previous study in Ref. |ll|] with ini- 
tial state |0, A/o,0) results in spin-mixing dynamics, due 
to which No was found to quickly reduce to some value 
without further oscillations or recovery. In our scheme, 
we find all exhibits collapse and revival patterns, so 
does Y as Y = N/3 — n - Even for the T 3 = block, we 
have seen redistribution of noise among the U-V com- 
ponents affects fluctuations in the isospin as well. The 




0.25 0.5 0.75 1 



X t 

a 

FIG. 3: Time-dependent squeezing parameter at (f> = 27r/3 
for N = 100 atoms, P = 1/3, and 6 = tt/2. 




FIG. 4: (a) The same as in Fig. 2 but for the optimized 
squeezing parameter; (b) the optimal angle <^ m i n which max- 
imizes squeezing as in (a). 

squeezing parameter 

^ {Ti)* + {T>)*' [ ' 

is analogous to ^ (^) but for isospin T' = R[4>]T after 
rotated around x-axis by an angle </>. Isospin squeezing is 
then characterized by £0 < 1. At <j) = 2it/3, this occurs 
after a very short time (see Fig. |^). It is especially 
interesting to note that £^ exhibits collapse and revival 
patterns. The optimal angle </> m i n for maximal squeezing 
(minimal Jll], |2^, ^] is shown in Fig. |J. It oscillates 
around its time-averaged value « 27r/3. In general, we 
find ^ achieves its minimum sooner and the minimum is 
smaller with decreasing values of 9 or increasing values 
of P . 

This effect is clearly unique to three-mode systems. 
In usual population spectroscopy (e.g. Ramsey type) 
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or in interferometry (e.g. of Mach-Zehnder type) for a 
two-mode system, particle partitioning noise and phase 
sensitivity can only be controlled by the modes involved 
directly. Here, the m/ = mode actually does not 
belong to the isospin group, yet it still influences the 
isospin noise properties. In contrast to the two-mode re- 
sult N± = J ± J z = N/2 ± 2T3, a three-mode system has 
N± = N/3 + (U 3 + V 3 ± 2T 3 )/2. A direct measurement 
of N + or N- will uncover all noise terms due to quan- 
tum correlations among the various spin components. A 
measurement of N + — N_, on the other hand, is similar 
to the two-mode case as the result is only affected by 
the noise in the isospin. When T3 = 0, the influence of 
the nif = mode population is reflected in the two-spin 
squeezing interaction between the U- and V-spins, which 
in turn also redistributes the noise in isospin. 

In Fig. |^, results of two-spin squeezing are shown for 
various initial Fock states \N- , No, N + ) of a spinor-1 con- 
densate. The lack of oscillations in Fig. 1(a) is due 
to non-oscillatory behavior of no for the particular ini- 
tial conditions used here. The solid curves are for the 
two-mode entanglement criterion of Ref . , valid only 
when No 3> N± . We see that when the initial states are 
such that N± modes are not near empty, the achievable 
two-spin or two mode squeezing essentially diminishes. 
However, there is a also turning point, when squeezing 
is again recovered if N± becomes significantly populated. 
Hence, we have found a new squeezing regime when the 
initial conditions are such that N± ^> Nq. The results 
are almost equivalent to the case No ^S> N± considered 
in Ref. [pT| . This new initial condition generates the 
two-mode entanglement via two spin squeezing between 
the U-V spin modes, i.e. between the holstein-Primakoff 
bosons. It should be noted that the two-mode entan- 
glement criterion in terms of spin squeezing parameters 
(£±) 2 has been derived for N > N± in Ref. 0. We 
show here that this criterion is also satisfied in the op- 
posite case of No <C N±. This observation emphasizes 
that the U-V squeezing criterion and the corresponding 
mode-entanglement can be sought for other initial condi- 
tions when the criterion of Ref. JTT[ | is no longer applica- 
ble. For that aim, we consider an initial state |25,0, 75) 
as shown in Fig. ^, where the U-V squeezing is indeed 
found. 



V. CONCLUSION 

We have provided a comprehensive treatment of quan- 
tum correlations in a spinor-1 condensate. Although 
no nonlinear interaction is apparent in the spinor con- 
densate Hamiltonian when single mode approximation 
is made, interesting quantum correlations do develop 
within subgroups of the SU(3) system. We have an- 
alyzed a spinor-1 condensate in terms of its T-, U-, 
and V-spin components. We have found and character- 
ized squeezing within one particular subgroup, similar 
to that of the isospin structure and we have numerically 



+ 



1200 




400 



FIG. 5: Time-dependent U-V squeezing parameter (dashed 
curve) and two-mode entanglement criterion (solid curve) for 
N = 100 atoms initially prepared in a Fock state of ?/<(0) = 
|iV_,iVo,iV+>: |0, 100,0) in (a), |1, 98, 1) in (b), |25, 50, 25} in 
(c), and |50,0,50) in (d). 



0.75 




!e ji/2 — 



7t/2 



0.002 0.004 0.006 0.008 0.01 

X 't 



FIG. 6: Same as Fig. ^ but now for the initial state ip(0) 
|25,0, 75). 



investigated its dynamics in terms of collapses and re- 
vivals. We have developed the U-V spin squeezing as a 
generalization of the often adopted spin (1/2) squeez- 
ing fill to two-spin squeezing. Its relation to mode- 
entanglement jllj in the Holstein-Primakoff representa- 
tion is also pointed out. We have presented new results 
for condensate fragmentation and spin-mixing phenom- 
ena in terms of the hypercharge symmetry and provided 
general phase-amplitude conditions for stationary states 
in the full quantum regime. 

In a typical experiment, a small magnetic field gradi- 
ent may be available |33|, which results in an effective 
Hamiltonian [§ H B = a(T + + T_) + /3T 3 2 -j B T 3 , instead 
of ([}]), with a, P, and 75 various renormalized parame- 
ters. In this case isospin squeezing still occurs through 
the one-axis twisting nonlinearity Jl3[ . 

Spin squeezing parameters can be measured directly by 
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the interferometry or Ramsey spectroscopy |L3]. Alter- 
natively, the isospin (T) squeezing in spinor-1 condensate 
can also be observed experimentally with light scattering. 
Using Raman coupled laser fields, an interaction of the 
type Hr = g(T + J- + h.c.) can be engineered (34|, [35[ , 
where J_ = y/2(aLa s + a} L aA) is an angular momentum 
operator, with aA,as,CLL the annihilation operators for 
anti-Stokes, Stokes, and pump photons. The interaction 
Hr allows for the mapping of spin correlations into pho- 
ton correlations as the total angular momentum T3 + J z 
is conserved. The solutions for J—(t) depend on the ini- 
tial conditions J_(0) and T_(0) ]34], |3f|. Therefore, the 
quadrature operators of scattered photons are directly re- 
lated to initial condensate spin quadratures and a homo- 



dyne measurement for Stokes parameters of the Raman 
field can reveal isospin squeezing []36f . 
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